(N 

o 

(N 
< 

in 



o 
o 



>: 

OS, 

vo, 
O' 
en; 

od 
O 

(N 



X 



Observation of fluctuation-induced tunneling conduction in micrometer-sized 
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Micrometer-sized Al/AlO^/Y tunnel junctions were fabricated by the electron-beam lithography technique. 
The thin (« 1.5-2 nm thickness) insulating AlO^ layer was grown on top of the Al base electrode by O2 glow 
discharge. The zero-bias conductances G(T) and the current-voltage characteristics of the junctions were 
measured in a wide temperature range 1.5-300 K. In addition to the direct tunneling conduction mechanism 
observed in low-G junctions, high-G junctions reveal a distinct charge transport process which manifests the 
thermally fluctuation-induced tunneling conduction (FITC) through short nanoconstrictions. We ascribe the 
experimental realization of the FITC mechanism to originating from the formations of "hot spots" (incomplete 
pinholes) in the AlO^ layer owing to large junction-barrier interfacial roughness. 

PACS numbers: 73.40.Rw, 73.40.Gk, 73.63.Rt 



I. INTRODUCTION 

Electron tunneling through a thin insulating layer (a 
potential barrier) in a metal-insulator-mctal (MIM) mul- 
tilaycred structure is one of the most fundamental re- 
search topics in condensed matter physics. It also lies at 
the heart of numerous solid-state devices, such as tun- 
nel diodesji Coulomb blockade thermometers^ Joseph- 
son junctions^ and memory elements based on magnetic 
tunnel junctions^ The functionality of a tunnel device 
relics heavily on the material properties of the intermedi- 
ate thin insulating layer. Usually, a weak insulating-like 
temperature dependence of the zero-bias junction con- 
ductance, G(T) = [dI(V,T)/dV] v ^o, as'described by 
the Simmons model^ is used to ascertain the quality and 
reliability of the insulating layer £ In the absence of any 
pinholes in the insulating layer [see the schematic dia- 
gram depicted in Fig. HJa)], the Simmons model predicts 
a quadratic temperature dependent G(T) oc T 2 law from 
low temperatures up to room temperature. 

In practice, a precise control of the material properties 
of the thin insulating layer during the junction fabrica- 
tion process is extremely difficult. Often, the junction 
yield varies sensitively with the oxide layer thickness i 7 - 
Furthermore, the thickness distribution of the oxide layer 
may be substantial even for those junctions grown under 
similar conditions^ As such, the G{T) behavior might 
differ substantially from the Simmons law in certain 
cases The key reason to explain these large varia- 
tions in the MIM quality can be ascribed to the common 
formation of junction-barrier interfacial roughnessj 11 ' 12 



a ' Electronic mail: yonghanlin@gmail.com 
b 'Elcctronic mail: jjlin@mail.nctu.cdu.tw 



Because the electron transmission probability increases 
exponentially with decreasing barrier width, the presence 
of any notable interfacial roughness thus can significantly 
affect the overall junction transport properties. 

In an MIM junction with marked interfacial roughness 
as schematically depicted in Fig. [T|b), the sharp points 
of the closest approaches (called "hot spots" hereafter—) 
between the two conducting regimes will no doubt play a 
predominant role in the electron tunneling. [We note that 
a hot spot is a precursor of a fully connected "pinhole," 
see Fig. [TJc) .] The current flowing through the hot spots 
between two large metal grains is theoretically found to 
depend sensitively not only on the applied bias voltage 
V a but also on the th ermal v oltage Vt across the hot 
spots, where Vt = ±\fk^T/C, with fee being the Boltz- 
mann constant, and C being the capacitance of the hot 
spotS i 14 i 15 (The plus/minus sign in Vt stresses the fact 
that Vt can be in the same/opposite direction to V a .) 
Owing to the smallncss of the hot spot, the magnitude 
of C is tiny and thus Vt is notable. The resulting effect 
of the total voltage Vtotai = V a + Vt gives rise to the so- 
called thermally fluctuation-induced tunneling conduc- 
tion (FITC) process through the junction. The FITC 
model predicts unique functional forms for the zero-bias 
G versus T dependence [Eq. (|2J)] as well as finite-bias 
current- voltage (I-V) characteristics [Eq. ((SJ] over a wide 
range of T from liquid-helium temperatures to room tem- 
perature. While this phenomenological FITC model has 
reasonably well described the G(T) behavior and I-V 
curves in a good number of experiments, such as nanowirc 
contacts ) 16 ' 17 the extracted barrier height is often too 
small (e.g., < few meV) and barrier width too large (e.g., 
> few tens nm). In order to resolve this puzzle, Xie 
and Sheng have recently retreated this problem by mi- 
croscopically calculating the electronic wave transmission 
through atomic-scale fine structures (short "nanocon- 
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FIG. 1. (a)-(c) Schematic diagrams depicting three types of 
insulating barrier: (a) a comparatively uniform oxide barrier, 
(b) an oxide barrier with two hot spots (close approaches) but 
no pinholes, and (c) an oxide barrier with two pinholes. The 
red dashed or solid curves in each case schematically depicts 
the electron transport process, (d) Schematic diagram show- 
ing the four- probe measurement method on a tunnel junction. 

strictions") connecting two large metallic grains J£ In 
particular, the relevant conduction channels they theo- 
retically considered are those with a transverse dimension 
slightly smaller than the half of the Fermi wavelength of 
the tunneling electron. By employing the Landauer for- 
mula, they obtained the G(T) attributes and I-V char- 
acteristics very similar to the original phenomenological 
FITC predictions. In other words, short nanoconstric- 
tions with a transverse dimension smaller than the cut- 
off wavelength (of a waveguide) can host thermally fluc- 
tuation induced tunneling conduction process at finite 
temperatures The extracted potential barrier height 
and width according to this microscopic model are pre- 
dicted to be much more realistic. 

On the experimental side, it is surprising that the pos- 
sible occurrence of the FITC-like mechanism in MIM 
multilaycred structures has never been reported, even 
though the formation of hot spots is very likely in certain 
junctions. In this paper, we present our observations of 
the FITC process governed G{T) features and I-V char- 
acteristics in micrometer-sized Al/AlO^/Y tunnel junc- 
tions in a wide temperature range 1.5-300 K. The ex- 
tracted parameters are examined in terms of the presence 
of hot spots in the AlO^ layer, which mimics the short 
nanoconstrictions in the microscopic FITC model. Our 
results arc presented below. 



II. EXPERIMENTAL METHOD 

The fabrication method of our planar tunnel junc- 
tions was described previously^ The designed geomet- 
rical junction area A<j (~ 20x20 /im 2 ) was defined by 



the electron-beam lithography technique. Both the bot- 
tom Al (« 30 nm thickness) and the top Y (« 100 nm 
thickness) electrodes were deposited by thermal evapora- 
tion. Prior to the deposition of the top Y electrode, an 
insulating AlO^ layer (« 1.5-2 nm thickness on average) 
was grown on the top surface of the Al base electrode 
by the O2 glow discharged Low angle x-ray diffraction 
studies revealed an amorphous structure of the oxide 
layer. The temperature dependent zero-bias resistance 
R(T) = 1/G(T) and the I-V curves of the junctions were 
measured by the four-probe method [Fig. Did)], using 
a Keithley K220 current source and a K182 voltmeter. 
A standard 4 He cryostat equipped with a calibrated Si 
diode thermometer was utilized for R(T) and I-V curve 
measurements between 1.5 and 300 K. The details of 
our sample configuration and measurement method had 
previously been discussed in Ref. [l9l 

We first discern the three possible types of junction- 
barrier interfacial roughness that can emerge in a fab- 
ricated MIM tunnel junction. Figure [l]a) depicts the 
normal case where the interfacial roughness is moder- 
ate. Electron tunneling occurs more or less evenly over 
the entire geometrical area A^, and the T dependence 
of the junction conductance G obeys the well-established 
Simmons law.— In the opposite case [Fig. [TJc)], pinholes 
may exist in the barrier, shorting the two electrodes, 
and cause overall metallic features of G(T). What is 
more interesting and will be the focus of this paper is 
the intermediate case [Fig. [lib)], where the interfacial 
roughness is notable but pinholes are just about to fully 
develop. Electron tunneling then occurs predominantly 
at the sharp points of the closest approaches of the top 
and bottom electrodes. Experimentally, the locations 
of such hot spots (incomplete pinholes) may be iden- 
tified by mapping the tunnel current over the lateral 
area A^, e.g., from conducting atomic force microscopy 
study at room temperature i 21 ' 22 Alternatively, the oc- 
currence of such hot spots can substantiate the FITC 
process in a wide T range from liquid-helium tempera- 
tures to room temperature ; 14 ' 15 as we demonstrate below 
in micrometer-sized Al/AlO^/Y tunnel junctions. 



III. RESULTS AND DISCUSSION 

According to the FITC model ) 14 ' 15 the small effective 
junction area A [i.e., the size of a hot spot depicted in 
Fig. Hlb), where A <C ^4d] can lead to large random ther- 
mal voltages Vt = ±^Jk^T/C that fluctuate across the 
narrow gap, which would in turn effectively lower and 
narrow the shape of the potential barrier. As a con- 
sequence, electron tunneling would be significantly in- 
fluenced, causing G(T) to increase exponentially with 
increasing T. (In the opposite limit of T —> 0, G be- 
comes constant, i.e., recovering the elastic tunneling be- 
havior, due to the gradually vanishing Vt with T.) Pre- 
viously, the FITC model has been applied to explain the 
G(T) behavior observed in a good number of conductor- 
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FIG. 2. (a) Zero-bias resistance R as a function of T and 
(b) the corresponding conductance G as a function of T 2 for 
the junctions A to D, as indicated. The straight lines drawn 
through junctions A and B in the top panel of (b) are least- 
squares fits to Eq. ([!]). 
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FIG. 3. (a) Log-log plot of R versus T for junctions C and D, 
as indicated. The solid curves are least-squares fits to Eq. 
The dashed curve through junction C is the theoretical pre- 
diction of Eq. (J2l) but is plotted by directly substituting the Ti 
and T values extracted from a(T) in Eq. ©. (b) Log(G/T) 
versus 1/T for junctions C and D between 25 and 300 K. Note 
that there exists no linear regime in any temperature regime. 



dielectric composite systems (granular films) 16 ' 17 ' 23 " — 
and nanowire contacts! 16 ' 17 Surprisingly, although the in- 
tcrfacial roughness must be very common in artificially 
fabricated MIM tunnel junctions ; 11 ' 12 ' 21 ' 22 ' 26 the possible 
manifestation of the FITC process through hot spots in 
the oxide layer has never been reported in the literature. 

In this work, a dozen of Al/A10 x /Y junctions have 
been fabricated and studied. In the following, we report 
the electrical-transport properties of four representative 
junctions. Two of them show the standard Simmons be- 
havior, while the other two manifest the FITC mecha- 
nism. Figure Ufa) displays the T dependence of R for the 
four junctions A to D, as indicated. In the two higher re- 
sistance junctions A and B, R increases only slightly with 
decreasing T and it progressively saturates at T < 100 K. 
In sharp contrast, in the two lower resistance junctions C 
and D, R increases exponentially with decreasing T. The 
relative resistance ratios [i?(1.5 K)-i?(300 K)]/i?(300 K) 
in the junctions C and D are larger than those in the 
junctions A and B. 

In the Simmons model, the zero-bias conductance, 
Gs{T), due to direct (elastic) tunneling of electrons 
through a potential barrier possesses a T 2 temperature 
dependence as given by^ 



G S (T) = G 



so 



T 



(1) 



where Gso and Tso are temperature independent param- 
eters. As shown in the top panel of Fig.^b), the G oc T 2 
law clearly holds for the junctions A and B in the wide 
temperature range 1.5-300 K. On the other hand, the 
bottom panel of Fig. Hfb) unambiguously indicates that 
the conductances of the junctions C and D significantly 
deviate from the Simmons law. 



In the phenomenological FITC theory, at small bias 
voltages (i.e., the ohmic I-V regime), the R(T) of a small 
junction is described b y 14 ' 15 



R(T) = Rcc exp 



T +T 



(2) 



where the parameter R^ depends only weakly on T, and 
the characteristic temperatures 



Ti = 



and 



16s 



7r(2m) 1 / 2 A:Be 2 u' 2 



(3) 



(4) 



where (f>o is the barrier height, w is the barrier width, 
So is the permittivity of vacuum, e r is the dielectric con- 
stant of the insulating barrier (for AI2O3; e r « 4.5-8.4, 
see Ref . l27t) . 2ttH is the Planck's constant, and m is the 
electronic mass. The characteristic thermal energy feeTi 
can be regarded as a measure of the energy required for 
an electron to cross the barrier, and Tq is the temperature 
well below which fluctuation effects become insignificant. 

Figure [3{a) shows a log-log plot of the measured R as 
a function of T for the junctions C and D. The symbols 
are the experimental data and the solid curves are the 
least-squares fits to Eq. J2]). The fitted Rao, Ti, and T 
values are listed in Table fl] The predictions of Eq. ^ are 
seen to well describe the data. (The dashed curve drawn 
through the junction C will be discussed below.) 

It might be conjectured that, if 4>o is small, the expo- 
nential increase of R with decreasing T in the junctions 
C and D could arise from the thermionic emission of elec- 
trons over the barrier i 2 ^ If this were the case, one would 
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TABLE I. Parameters for junctions C and D. The values of 
w and 4>o were extracted by using the geometrical area Aa (~ 
20 x 20 ^im 2 ), and the dielectric constant e r (Al203) ps 6. 



Junction 


Method 


Roo 


Ti 


To 


w 


4>o 






(kfi) 


(K) 


(K) 


M 


(meV) 


C 


R(T) 


7.66 


815 


380 


37.7 


0.050 




/h(V h ) 


6.46 


915 


391 


39.5 


0.054 


D 


R(T) 


11.4 


21.5 


46.5 


22.9 


0.006 



expect the total conductance to be given by the sum of 
the Simmons conductance and the thermionic-emission 
conductance G t (T), i.e., G(T) = G S (T) + G t (T), where 
G t (T) = G to Texp(-T t o/T), with G t0 and T t0 being 
two temperature independent parameters. At low tem- 
peratures, Gs(T) should dominate over Gt(T). As T 
increases, Gt(T) would become progressively important 
and might eventually dominate over Gs(T). Figure EJb) 
shows a plot of log(G/T) versus 1/T for the junctions C 
and D. Clearly, there does not exist any linear regime be- 
tween 25 and 300 K. Therefore, the thermionic-emission 
conduction process can not account for our measured G 
versus T dependence. 

We return to the FITC mechanism. As recently dis- 
cussed in Ref. [l7| by two of the authors, one may alterna- 
tively extract the values of T\ and T by measuring the 
non-ohmic I-V characteristics of a given junction at high 
bias voltages Vh- Under such conditions, the phenomeno- 
logical FITC theory predicts a nonlinear I-V curve given 
by^ 



h(Vh,T) = 7 hs exp 



(5) 

where the saturation current i^s and the critical voltage 
Vhc arc two parameters which depend only weakly on T. 
The parameter a(T) = T\/(Tq + T) describes the tem- 
perature effect on the I-V curves. The value of a(T) at 
each temperature can be extracted by fitting the mea- 
sured Jh(Vh) curve to Eq. ([5]). Then, the values of Ti 
and To can be inferred by fitting a(T) to the expression 
T\/{T Q + T). Note that with the two sets of Ti and 
To values deduced from the two complementary manners 
[Eq. ([2} and Eq. ((3J)] , one may perform a quantitative self- 
consistency check of the FITC theoretical predictions. 

Figures |4ja) shows the measured nonlinear I-V curves 
for the junction C at four T values. The symbols are 
the experimental data and the solid curves are the least- 
squares fits to Eq. ([5]). A good agreement between the 
theory and experiment is evident. Our extracted values 
of a as a function of T are plotted in Fig. 0Jb), together 
with the least-squares fits (solid curves) to the expression 
Ti/(T +T). The Ti and T values thus inferred are listed 
in Table |U The insets of Fig SJb) plot the variations of 
our fitted Jh s and 14, c values with T. Note that they are 
essentially temperature independent, as predicted by the 
FITC theory. 




(b) 




4 8 12 16 20 

K(mV) 



FIG. 4. Current- voltage characteristics of junction C: (a) non- 
linear I-V curves at four T values. The solid curves are least- 
squares fits to Eq. ([5]). For clarity, the data for 15, 30, and 
45 K have been shifted up by multiplying factors of 1.2, 1.4, 
and 1.6, respectively. Inset: the I-V curve at 1.5 K in double- 
linear scales, (b) The parameter a in Eq. ((SJ as a function 
of T. The solid curve is the least-squares fits to Ti/(To +T). 
The dashed curve is the theoretical prediction of Ti / (To + T) 
but is plotted by directly substituting the Ti and To values 
extracted from the R(T) fit to Eq. Insets: (top) variation 
of Ihs and (bottom) variation of Vhc with T. 



We now substitute those T\ and T values extracted 
from the R{T) fit [Eq. ©] into the expression Ti/(T +T) 
to compute a(T) for the junction C. The calculated re- 
sult (dashed curve) is plotted in Fig. EJb), which is seen 
to describe the a(T) behavior reasonably well. In turn, in 
Fig. [3fa) we plot the prediction of Eq. ([2]) by directly sub- 
stituting those Ti and To values inferred from the a(T) fit 
through Eq. (while allowing R^ as the sole adjustable 
parameter, whose value is listed in Table [J). The result 
(dashed curve) is clearly in excellent agreement with the 
experimental data. Numerically, our T\ and To values 
extracted from the two methods differ by an amount of 
ps 10%, which is very satisfactory and notably smaller 
than those (ps 30-40%) reported in previous studies. 17 ' 24 
This close consistency provides a strong support for the 
pivotal relevance of the FITC mechanism occurring in 
our micrometer-sized tunnel junctions. 

We would like to comment on our inferred values of 
the barrier parameters </>o and w. Inspection of Table U 
indicates that the </>o values are relatively small (< 0.05 
meV), while the w values are relatively large (> 20 nm). 
These seemingly nonphysical values can be reconciled by 
taking into account the possible formation of hot spots 
as those depicted in Fig. [ljb). Since the hot-spot area 
A -C , our using the geometrical area A^ to extract 4>o 
and w from the inferred T\ and To values through Eqs. ([3]) 
and ([4| hence led to huge underestimate (overestimate) 
of 4>q (w). Although the actual area A of a hot spot can 
not be precisely known for a given tunnel junction, it has 
previously been pointed out that the typical size of a hot 
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spot can be about a few nanometers. 21 ' 22 Then, taking 
Aw 2x2 nm 2 , one obtains values of </>o w 80 (10) meV 
and w « 0.9 (0.6) nm for the junction C (D). These 
values are much acceptable. 

In fact, the relevant hot-spot area could even be re- 
duced to the sub-nanometer scale. In order to search 
for a microscopic insight and to obtain meaningful pa- 
rameter values, Xie and Shengi^ have recently consid- 
ered the electron tunneling through hnite segment (s) of 
nanoconstrictions whose transverse dimension is less than 
the half of the electronic Fermi wavelength Af- In this 
new theoretical proposal, the insulating gap in the origi- 
nal FITC models is replaced by short and very narrow 
constrictions (which behave like waveguides with trans- 
verse dimensions slightly shorter than the cutoff wave- 
length Af/2). By applying the Landauer formula and 
taking into account the effects of thermally induced volt- 
age fluctuations across the nanoconstrictions, Xie and 
Sheng obtained a G versus T dependence very similar to 
the original FITC prediction [Eq. ©] but with more re- 
alistic parameter values. Following this new theoretical 
idea and taking Am 0.2x0.2 nm 2 [note that Af(A1) « 
0.36 nm], we obtain </>o « 0.5 eV and w ~ 0.4 nm for the 
junction C. These values are fairly realistic. Therefore, 
the manifestation of the FITC mechanism in micrometer- 
sized tunnel junctions implies the existence of such fine 
structures (hot spots or short nanoconstrictions) in the 
oxide layer. Electron tunneling through these fine struc- 
tures governs the overall G(T) behavior. [There could 
be a few such fine structures coexisting in a given MIM 
junction. However, the measured G(T) might be pre- 
dominated by only one of them, owing to the exponen- 
tial dependence of the electron transmissivity on the fine- 
structure geometries.] 

In a series of ballistic electron emission microscopy 
(BEEM) and scanning tunneling spectroscopy (STS) 
studies, Buhrman and coworkers^^ have systematically 
shown that in an ultrathin (< 1 nm) and not fully ox- 
idized A10 x layer (which was not covered with a top 
electrode), there could possibly exist low-energy elec- 
tron channels that could provide low-voltage "leakage" 
through the barrier. In the present study, our AlO^ lay- 
ers are relatively thicker (1.5-2 nm) and fully oxidized by 
the O2 glow discharge, as mentioned in Sec. II. There- 
fore, we think that the features of the FITC conduction 
process observed in our Al/A10 x /Y junctions are not as- 
sociated with the possible presence of low-energy electron 
channels that extended through the thin disordered oxide 
layer. 

Theoretically, the crossover from a tunneling regime 
to an incoherent thermal-activation regime in a cor- 
related barrier has been investigated by Freericks and 
coauthors.— They have shown gcnerically that the 
crossover takes place around a characteristic temperature 
which approximately equals a generalized Thouless en- 
ergy for the barrier. Under certain parameter values, the 
predicted temperature behavior of the junction resistance 
(or, the specific resistance RA^) looks nominally similar 



to the temperature dependence shown in Fig. |3](a). Un- 
fortunately, the theory of Freericks does not provide a 
functional form for comparison with experiment. This 
issue deserves further studies. 



IV. CONCLUSION 

We have observed the fluctuation-induced tunneling 
conduction mechanism in micromctcr-sized Al/AlO^/Y 
tunnel junctions in a wide temperature range 1.5-300 K. 
The manifestation of the FITC process at the microm- 
eter scale reflects the existence of large junction-barrier 
interfacial roughness in the thin oxide layer, where elec- 
tron tunneling occurs at the sharp points of the closest 
approaches of the top and bottom electrodes. The forma- 
tion of a few hot spots, but not pinholes, predominantly 
governs the conductance versus temperature behavior as 
well as the current-voltage characteristics. Our results 
can have important bearing on the reliability and func- 
tionality of solid-state tunnel devices. 
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